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We present generic scaling laws relating spreading criti- 
cal exponents and avalanche exponents (in the sense of self- 
organized criticality) in general systems with absorbing states. 
Using these scaling laws we present a collection of the state- 
of-the-art exponents for directed percolation, dynamical per- 
colation and other universality classes. This collection of re- 
sults should help to elucidate the connections of self-organized 
criticality and systems with absorbing states. In particular, 
some non- universality in avalanche exponents is predicted for 
systems with many absorbing states. 

PACS numbers: 64.60.Lx, 05.40. +j, 05.70.Ln 



Directed percolation (DP) is broadly recognized as the 
paradigmatic example of systems exhibiting a transition 
from an active to an absorbing phase DP criti- 

cal behavior appears in a vast array of systems, among 
others chemical reaction-diffusion models of catalysis || , 
the contact process damage spreading transitions 

||, pinning of driven interfaces in random media 
roughening transitions in one-dimensional systems 
and Reggeon field theory || . This universality class has 
proven very robust with respect to the introduction of 
microscopic changes, and many apparently different sys- 
tems share the same critical "epidemic" or "spreading" 
H and "bulk" exponents Nevertheless, examples 

of a system exhibiting a transition to an absorbing state 
outside the DP class have been identified in recent years. 
Some examples are: 

• Systems with two symmetric absorbing states or, 
what is equivalent in many cases, systems in which 
the parity of the number of particles is conserved 



• Systems with an infinite number of absorbing 
states, which exhibit nonuniversal spreading expo- 
nents HQ. 

• Systems in which the dynamics is limited to the 
interface between active and absorbing regions. 
These are in the class of the exactly solvable voter 
model and compact directed percolation Jl5| . 

• Some models of epidemics with immunization (no 
reinfection) |L6]]. These belong to the so-called dy- 
namic percolation class; the final set of immune 
sites at criticality is a percolation cluster. 



Recently, connections between self-organized criticality 
(SOC) and systems with absorbing states have attracted 
much attention. For example, there has been a debate 
on whether the extremal Bak-Sneppen for punctuated 
evolution jlT]] and certain variants are related to DP [|f8] . 
It has also been argued that sandpile models [^9| share a 
number of features with systems having many absorbing 
states p0| |) and certain self-organized forest-fire models 
are related to dynamical percolation 21 1. 

In self-organized models the so called avalanche expo- 
nents are customarily determined. Surprisingly, in spite 
of their obvious similarities, the general connections be- 
tween spreading and avalanche exponents have not, to 
the best of our knowledge, been given explicitly for gen- 
eral systems with absorbing states. Establishing the gen- 
eral scaling laws relating avalanche and spreading expo- 
nents in systems with absorbing states is the main goal 
of what follows. This will allow us to put together many 
different scaling relations and exponent values, presently 
quite dispersed, and sometimes difficult to find in the 
literature, and should facilitate progress in this field. 

Let us first define in detail spreading and avalanche 
critical exponents. The most accurate determination 
of the critical point of systems with absorbing states 
comes from "epidemic" or "spreading" experiments 
In these, a small perturbation (localized activity) is cre- 
ated at the origin of an otherwise absorbing configura- 
tion, leading to a spread of activity. In spreading exper- 
iments, it is customary to measure the number of par- 
ticles, averaged over all runs (including those that have 
reached the absorbing state) N(t), the survival probabil- 
ity P(t), and the mean-squared deviation from the origin 
R 2 (t) . At criticality these magnitudes scale as 



N(t) ~ e>; p{t) ~ r 



R 2 {t)~t z 



(1) 



where r/, 5 and z p2fl are the spreading exponents. 

Once the critical point has been located accurately 
all the remaining standard critical indices can be esti- 
mated. For reference we show in Table I a compilation of 
the state-of-the-art values for the usual exponents in di- 
rected percolation, corresponding to a synthesis of the 
best series expansion and simulation results. For the 
sake of completeness let us give here their respective def- 
initions. Calling A the distance to the critical point 
in terms of the reduced control parameter, p the order 
parameter, £j_ (£||) the characteristic length (time), h 
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an external field coupled to p, and \ = L d vai(p) the 
static "susceptibility", we have: p ~ A' 5 , p ~ /i 1 /" 5 ' 1 at 
the critical point, x ~ A" 7 ', f ± - A - "- 1 -, C|| ~ A~' y ii , 
Poo = linif^oo P(t) ~ A' 3 , and ~ at the critical 
point. 

^From the whole set of exponents that can be defined 
in DP, only three are independent; the rest can be deter- 
mined using well-known scaling relations (see appendix). 
In certain systems possessing an infinite number of ab- 
sorbing states a fourth independent critical ex- 
ponent has to be introduced [^3 24 1. This is due to the 
fact that the exponent 6 (which in DP coincides with 9) is 
non-universal and depends on the nature of the absorbing 
state in which the epidemic spreads [ ^3|j24| . (Similarly, 
the exponent /?', normally identical to /3, varies along 
with 5 in such systems [p3|.) 

On the other hand, studies of avalanche transport em- 
ploy a different definition of the spread of activity. For in- 
stance, in the prototypical sandpile model Jig) ] avalanches 
are obtained by adding one sand grain to a stable or 
absorbing configuration. In this way the system jumps 
among absorbing configurations via avalanche-like rear- 
rangements. The following quantities and associated ex- 
ponents are usually measured: 



P(s) ~ s- T f(s/s c ) 



Sc 

(s) 



(2) 
(3) 
(4) 



where s is the size of an avalanche, (total number of ac- 
tive or toppling sites), P(s) the associated probability 
distribution, s c the cut-off size, (s) the mean size, and e 
represents the temperature-like variable associated with 
the process: A for contact process or DP, the dissipa- 
tion rate for sandpiles, F — F c in driven-interface models 
(here F is the driving force). If e = 0, the characteristic 
length is defined by the system size L through the scaling 
relation s c ~ L D . Analogously, the following exponents 
associated with the duration t are also measured 



P(t) 

to 

(*} 



t- T 'g{t/t c ) 
f -i/ff< 



"7t 



(5) 
(6) 
(7) 



Let us now provide the general scaling laws relating 
avalanche and spreading exponents in systems with ab- 
sorbing states. 

^From the definitions of r\ and 5, it is evident that 
the total number of particles in surviving runs goes like 
N s ~ t n+s , and therefore its time integral is governed by 
the exponent 1 + 77 + 5. Thus an avalanche that dies at 
time t has a typical size s ~ 

f i+*7+<5. The probability to 
die between times t and t + dt scales as D(t) ~ t~ d ~ 1 dt. 
Observe that the time is defined in such a way that after 
a 'toppling' (updating of a given site), it is increased by 
At = 1/N s (t) [[25). Therefore the number of topplings 
per unit time is N s (t). To express r, a and 7 as func- 
tions of the spreading exponents, let us consider a specific 



avalanche size, say si. An avalanche of size si can have 
different durations, since t and s are not related in a de- 
terministic way, i.e., 



P(si) 



dtP(sx\t)D(t) 



(8) 



where t\ and £2 are the minimum and maximum times 
compatible with si, and P(s\t) is the conditional proba- 
bility of an avalanche having size s, given it dies at time t. 
P(s\t) is bell-shaped, with its maximum at t ~ s 1 /( 1 + , '+' 5 ) 
M. Writing P(s\t) = £-( 1 +'?+' 5 ) F(s/t 1+ ^ +s ), where 
F[u), the (normalized) scaling function, is nonsingular, 
we have, on changing variables 

P(s) = s -(i+v+M)/(i+v+s) J duu mi+v+s) F ^ 5 (9) 



in other words, 



l + rj + d ' 

With t c ~ e^n , and using s c ~ t 1 c +n+6 , we have 
l/a = v\\ (l + ?7 + 5). 



(10) 



fll) 



The remaining exponent 7 and the fractal dimension 
D can be determined using the relations 7 = (2 — t)/<7 = 
v \\ (1 + 7 ?) (this last equation for 7 has already been found 
by other authors ^7j), and the standard relation D = 
l/(av±). Following a very similar derivation to the one 
just presented, one can easily determine also the following 
scaling relations for the exponents associated with P(t): 



n = l + 5 

o~ t = f(l + V + 5) = 1/V[| 
Jt =(2-T t )/a t = p»(l-S). 



(12) 



Let us derive explicitly the expression for j t . In DP, 
we have, for A < 0, or for A = and finite L, the 
scaling form for the survival probability P(t) ~ t~ s e~ t l tc 
with t c ~ |A| _1/ ii, or t c - L v \^ u± . Since the probability 
density for dying at time t is —dP(t)/dt, we can write 



(t) = I t-[-P(t)]dt 



r 



t 



'I. 

1-5 



u- 5 e~ u du. 



(13) 



So (t) ~ lAI^ 1 -^!! , giving the value of j t . In the finite- 
size case, (t) ~ jj-^^w/"^ (observe that to and wo are 
unimportant lower cut-offs). 

All the scaling relations derived so far are general, and 
valid for all systems with absorbing states. Specific scal- 
ing relations for systems in the DP class can be written 
using the well known relation || 77 + 25 = dz/2. Using 
the best values for the spreading exponents in DP taken 
from the bibliography (Table I), we obtain the values of 
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the avalanche exponents for DP in different dimensions 
(they are also summarized in table I). 

Applying our general relations to other classes of mod- 
els, we obtain: 

• For models with parity conservation, using the 
known result for spreading exponents |llf| , we pre- 
dict t w 1.22, a w 0.24, 7 » 3.25, r t « 1.28, 
at ~ 0.31 and 7< w 2.33 in d = 1 [£8| an d mean 
field values above that dimension. These results 
have also been derived and numerically tested in 

• In systems with many absorbing states, a general- 
ized hyperscaling relation has to be introduced, due 
to the fact that in this case the exponents 5 and rj 
are non-universal and therefore 5 =/= 9 in general 
(on the other hand the combination r\ + 8 retains 
its DP value}. The generalized scaling law for these 
systems is f] + 6 + 9 = dz/2. Applying our 
scaling laws we predict non-universal values of t, 7, 
Tt and "ft for systems with many absorbing states; 
i.e. if experiments are performed on a fixed envi- 
ronment, the results depend upon the environment 
itself. Recently this kind of non-universality has 
been observed in the class of sandpile models with 
fixed energy Q . 

• For models in the CDP class we have r = 4/3, 
a = 2/3, 7 = 1, T t = 3/2, a t = 1 and 7t = 1/2 in 
d = 1 and mean field values in d = 2 and above. 

• For dynamical percolation we can take advantage of 
our scaling laws, using them the other way around, 
i.e., using the well known avalanche (cluster) expo- 
nents for standard percolation ]3l|,|37| permits us to 
determine the spreading exponents [p4[ with good 
accuracy. In table II, we report a collection of ex- 
ponent values in d = 2, 3 and 6 spatial dimensions. 

In summary, we have presented the general scaling re- 
lations that rule general systems with absorbing states, 
and present a collection of exponent values that can be 
useful as a reference. We believe that this coherent 
derivation and collection of otherwise scattered scaling 
laws and exponent values may facilitate progress in draw- 
ing connections and similarities in many systems which 
show absorbing states and avalanche behavior. 
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After completion of this work we became aware of a 
recent paper by Lauritsen et al. |38|, in which very similar 



they find in that case S ^ 9. A direct consequence is that, 
as discussed here, some avalanche exponents do not take 
their corresponding DP values. Some other interesting 
scaling relations can be found in J3^] . 

APPENDIX: Scaling relations for DP 

Here we present a collection of scaling laws for the DP 
universality class |j]||[j)< 



rj + 5 + 9 = dz/2 



S h 
5 

i 

(3- 

z ■■ 

D : 

7 : 



7 — ^ = dv x — 2/3 

Bv \\ 

8v \\ 

2v±./v\\ 

l/(av_ L )=d+(v\\ -0)/ux 
(2 - T)ja = dv± + - 2/3. 



(14) 



Observe that not all of these relations are independent. 



scaling relations to the ones proposed here are derived for 
directed percolation in the presence of an absorbing wall. 
In particular, as in systems with many absorbing states 
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Critical exponents for Directed Percolation. Exponents cal- 
culated by using scaling relations contained in this paper are 
reported in the lower part. ° See |p7}| ; b (32); c obtained using 
(3 = 6v u ; d ||; e obtained using v±_ = zv ]} /2; f §|; 9 ||; 
h Where not reported uncertainties are in the last digit. 
For d — 4 we report the exact mean field values. 
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Critical exponents for dynamical percolation. Exponents cal- 
culated by using scaling relations contained in this paper are 
reported in the lower part. The rest of exponents values are 
from P?| . Where not reported uncertainties are in the last 
digit. For d = 2, values expressed as fractions refer to exact 
results . For d = 6 we report the exact mean field values. 
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